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QJ , Abstract 



We study the Chern-Simons number diffusion rate in the (l-l-l)-diniensional 
lattice Abehan Higgs model at temperatures much higher than, as well as com- 
;J] ' parable to, the sphaleron energy. It is found that in the high-temperature limit 

^ the rate is likely to grow as power 2/3 of the temperature. In the intermediate- 

temperature regime, our numerical simulations show that very weak temperature 
dependence of the rate, found in previous work, persists at smaller lattice spac- 
ings. We discuss possibilities of relating the observed behavior of the rate to static 
finite-temperature properties of the model. 



*Address after 1 September 1994: Niels Bohr Institute, Blegdamsvej 17, DK-2100 Copenhagen, 
Denmark. 



The Abelian Higgs model (AHM) in 1+1 dimensions is one of the simplest yet non- 
trivial field-theoretic analogs of the Standard Model (SM) as far as fermion number 
nonconservation at high temperature is concerned. Clearly, there is little hope of quan- 
titatively understanding this phenomenon in SM if a similar goal cannot be achieved 
for AHM, which lends itself much more easily to both analytical and numerical studies. 
This explains the considerable amount of work done on this aspect of AHM in recent 
years ]I|,g,|,|,|, 1,0. 

Prior to the present work, the status of knowledge on thermally induced fermion- 
number violation in AHM has been as follows. Much like in SM, the dynamics of the 
Chern-Simons number A^cs (and thus of the fermion number) in AHM can reliably be 
studied in the classical approximation at temperatures well above the gauge- (mw) and 
Higgs- [rriH) boson mass scales. The relative simplicity of AHM allows to give an accurate 
estimate of the A^cs diffusion rate at temperatures T low compared to the energy Egph of 
a sphaleron saddle-point configuration separating vacua labeled by consecutive integer 
values of A^cs- Not surprisingly, the rate F exhibits thermal-activation behavior 

7 

r = im'/' (^) ' exp (-^) , m^, mn « T « E,,^, (1) 

where ^ = 2my^/m\. Unlike the SM case, in AHM both the sphaleron energy and 
the / prefactor are known analytically. Numerical real-time simulations of AHM at low 
temperatures, including those of the current work, indeed agree with (|Tp |I|, ^, |^. 

Surprising as it may sound, at the opposite end of thermal scale, T ^ -Egph, less was 
known about the rate in AHM than in SM. In SM, the only dimensionful quantity relevant 
at such high temperatures is the temperature itself, hence, on dimensional grounds, the 
rate is expected to behave as T^. In AHM, however, both the gauge and the quartic 
Higss coupling constants are dimensionful, and therefore, while a power-law behavior was 
generally believed to be satisfied P, § , the corresponding power could not be established 
from dimensional considerations. Nevertheless, a simple scaling argument given in the 
following indicates that this power is 2/3. This power law is also consistent with our 
numerical results. 

As usually is the case, the intermediate-temperature regime (T comparable to -Egph) 
is the least accessible to analytical estimates. Numerical work 0, |] revealed a surprising 
feature: at temperatures 5 to 7 times less than E^ph the rate practically stops grow- 
ing and remains roughly constant until the sphaleron scale of temperatures is reached. 
Our simulations show that this phenomenon persists at smaller lattice spacings and is 
therefore not a lattice artifact. 

We shall now elaborate on the points just made, discussing first the physics and 
afterwards explaining our numerical methods. We begin by giving a scaling argument 
in favor of the T^/^ behavior of the rate at high temperatures. Consider the lattice- 
regularized classical AHM whose temporal-gauge Hamiltonian can, with a suitable choice 



of units, be written in a dimensionless form p, ^ : 

H{v, o) = ^ E[e4.+i + K.f + y (I'^.f - ^y + I'^.+i - e^'^^--0,r], (2) 

where ttj and -Ej,j+i are canonically conjugate momenta of the complex scalar (Higgs) 
field (pj and the spatial component Ajj+i of the gauge field, residing respectively on sites 
and links of an A^-site chain. The lattice spacing is a. Periodic boundary conditions 
are assumed. The reader can verify by direct substitution that if 0j(t), Tfjit), Ajj^i(t), 
Ejj^i{t) are solutions at time t of the equations of motion following from H{v,a) with 
H{v,a) = £, then c(j)j{ct), C7ij{ct), cAjj+i{ct), cEjj^i{ct) play the same role at time 
t for H{cv, a/c) with H{cv, a/c) = c'E. Because of this property, scahng laws can be 
derived for both static and dynamical observables in AHM. Consider in particular the 
Chern-Simons number defined as 

^cs(t) = |^EA,,+i(t). (3) 



Its diffusion rate is 



V{a,v,T) = lim^,,_oo^((iVcs(t) - iVcs(0))')T. (4) 



By ()t we mean integration over all the initial values of dynamical variables [10| with 
the Boltzmann weight 

exp(-^^)n/(^.H). (5) 

where the delta-function factors impose a set of Gauss' constraints 

C,{a) ^ -(E,-,+i - E,_i,,) - Im(7r,0,) = 0. (6) 

Multiplying all the integration variables {i.e., all the initial values of fields and momenta) 
in (I) by c and using the scaling property of H{v, a) under this transformation we obtain 
a simple scaling rule for the rate: 

r (-, cv, c^t\ = c^T{a, V, T). (7) 

Differentiating with respect to c at c = 1 gives a differential equation for F as a function 
of a,v,T: 

{vdy + 3TdT -ada-2)T = 0, (8) 

whose most general solution is 

r = Tig Ut, Pi , (9) 



0^025 



0.020 



O 

^~ 

CO 

m 



0.015 



0.010 




0.02 



0.03 



0.04 



0.05 



0.06 



0.07 



aT^' 



Figure 1: g{a^T, 0) measured at small values of aT^'^. 



where g{x, y) is an arbitrary function of its two arguments. In particular, since E'sph oc w^, 
(0), derived for the continuum theory, shows that g{0,y) oc y'^^^ exp{—y) for y ^^ 1. As 
for the temperatures far above the sphaleron scale, the continuum limit of the rate, if it 
exists, is ^'(0, 0)T^/^ as claimed. Both result (|I|) and numerical evidence strongly suggest 
that the rate does indeed have a finite continuum limit for any nonzero v^/T. It is also 
plausible that the rate becomes insensitive to the details of the Higgs potential at energy 
scales far below T, i.e., T becomes independent of v ii v^ <^ T. We therefore expect 
< ^(0,0) < oo. 

We verified the scaling rule (|^ by a numerical simulation. Since in this part of the 
work we are primarily interested in the high-temperature regime, we simplify matters by 
setting V = 0. The rate measurement then gives 

T{a = 1/2, T = 1 



r(a = l,T = l/8) 



4. ±0.4, 



(10) 



in agreement with (|^). Next, we performed a series of rate measurements closer to the 
continuum limit of the v = theory in order to estimate g{0,0). The measurements, 
done at lattice spacings 1/16 and 1/8, span a range of values of the scaling-invariant 
parameter aT^f^ between 0.028 and 0.065. The results, presented in Figure |I|, indicate 
that g{x, 0) approaches a finite value as x — >■ 0. 

In this work we are concerned with classical AHM. Note, however, that if the classical 
approximation for the rate is valid in quantum AHM at very high temperatures, so is 
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Figure 2: Measured rate dependence on the inverse temperature (3. The sohd hne is 
the activation-theory prediction. The dashed hne shows the T^/^ law for the high- 
temperature theory, with 5'(0,0) estimated at /5 = 11, a = 1/16. 



the T^/^ law. Indeed, for consistency of the classical approximation the classical lattice 
spacing is chosen inversely proportional to the temperature |jl|, |^ , yielding a^T -^ 
required in order for the T^/'^ law to hold. 



as 



Next, we studied the intermediate-temperature properties of the rate. As in previous 
work on the subject, we fixed f = 1, resulting in i^sph ~ 0.94. The inverse-temperature 
range of interest is then 1 < /5 < 5, where an extremely weak temperature dependence of 
the rate was found in previous work 0, Q . By simulating the theory closer to continuum 
limit (a = 1/16 compared to a > 1/8 in the earlier work), we verified that this feature 
is not a lattice artifact: as Figure H shows, the plateau region becomes wider and more 
pronounced as the lattice spacing shrinks. At still lower temperatures, /3 > 7, the 
lattice-spacing dependence of the rate rapidly diminishes. We find a good agreement 
between our low-temperature rate measurements and those of |Q, performed at larger 
lattice spacings. Both earlier and current results approach the asymptotic form (Q) with 
decreasing temperature. 

Let us summarize the properties of Nqs diffusion rate in AHM as they emerge from 
this and previous works. At low (T -C £^sph) the rate follows the activation-theory predic- 
tion (p and is exponentially suppressed by a Boltzmann factor involving the sphaleron 
barrier height. At extremely high temperatures (T ^ -Esph) the rate, if it has a finite con- 
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Figure 3: Relative variation of the density of states with energy per degree of freedom 
in AHM (sohd hne) and in a free- field theory (dashed line). 



tinuum limit, should grow as T^/^. Our simulations do point to a finite continuum limit 
of the high-temperature rate. The crossover between the two extreme regimes involves a 
temperature interval characterized by a very slow or null rate growth. A complete theory 
of A^cs diffusion should contain a derivation of the T"^^^ law at high temperatures and a 
determination of g{0, 0) from first principles. It also should explain how corrections to 
either high- or low-temperature rate lead to the observed interpolation between the two 
extremes. In the absence of such a theory the best we can do is look for clues that might 
help construct one. 

To this end, we examined a possibility that the reasons for the rate suppression 
compared to (|l]) at intermediate temperatures are entropic in nature |^. Namely, the 
derivation of (|l]) is based on approximating the AHM phase space by that of linear 
excitations about the ground state. Starting at very low temperatures T, where it is 
definitely valid, this approximation will only hold with increasing T as long as the free- 
field relation 

5^1og(v^mp(^(T)))=5^ = i (U) 

approximately holds. Here E{T) is the average energy of the system at temperature T, 
e(T) is the average energy per degree of freedom, and p{E) is the corresponding density 
of states. We determined the left-hand side of (0) using the binning method fill. As 
Figure |^ shows, the variation of the AHM density of states hardly differs from, and 
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Figure 4: Temperature dependence of (|0P)t in t; = 1 (squares) and v = (crosses) 
theories. 



even appears somewhat smaller than, that of a free-field theory, even as the energy per 
degree of freedom (approximately equal to T/2 throughout the range of temperatures 
in question) approaches the sphaleron scale. We therefore conclude that the observed 
slowdown of the rate growth is unlikely to be entropy-driven. 

In search for other clues we observed that in the intermediate range of temperatures 
(a) the average squared expectation value (0^)t of the scalar field is close to a minimum 
(Figures and ^); (b) the rate in the theory with f ^ = is considerably below that of 
f ^ = 1 theory, while (c) (|0P)t is considerably greater in the f ^ = 1 theory than in the 
t>^ = one. Put together, these pieces of empirical evidence suggest that the rate might 
be closely related to (|0P)t- Moreover, repeating our scaling argument for (|0p)r, we 
obtain an analog of (^: 

3\ 
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a^T,^ 
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where, again, we expect < h{0, 0) < oo. As our measurements (Table |^) show, 
h{a^T,0) = {\(J)\'^)t/tI varies by only about 5% for a^T between 1.9 x 10"^ and 1.25 x 



10~^. We conclude that, in all likelihood, (|0P)t oc T^^^ in the continuum limit of 
high-temperature AHM. As a consequence, we also expect F oc (|0P)t in this limit. It 
therefore appears worthwhile to explore the possibility that at intermediate and high 
temperatures the time scale for A^cs diffusion is primarily set by the value of (|0P)t- 



a^T X 10^ 


12.6 


8.20 


4.93 


3.52 


2.73 


2.24 


1.89 


(10 ^)t 


0.615 
±0.005 


0.605 
±0.004 


0.600 
±0.004 


0.594 
±0.004 


0.595 
±0.007 


0.599 
±0.004 


0.588 
±0.004 



Table 1: Average squared magnitude of the scalar field as a function of o?T in the f = 
theory. 

This seems plausible given the equation of motion for the Chern-Simons number written 
in the form 

Ncs = ^ 51 PjPj+^ sin 6jj+i = — (|0Ht ^ sin fo^j+i ± corrections, (13) 

where aj = arg0j, 0j = pj exp(iaj), and bj = aj+i — otj — aAj,j + 1. 

We conclude by explaining simulation methods used in this work. We chose to deter- 
mine r as given by @ by generating a thermal ensemble of initial conditions for lattice 
AHM. We then let each of these initial field configurations evolve according to the AHM 
Hamiltonian equations of motion. The rate was obtained by averaging F of (^) both 
over the ensemble of initial conditions and within individual Hamiltonian trajectories of 
the system ||^. Initial configurations were generated by applying standard Monte-Carlo 
techniques to the theory written in terms of gauge-invariant variables |^ : 

j 

where pj and 6jj+i are defined as in (|13|) , ejj+i = Ejj+i/a, and tTj is the canonical mo- 
mentum of Pj. More specifically, we used the Metropolis algorithm to update the p and 
b variables and the heat-bath algorithm to update e. The fourth canonical variable, 7Tp, 
is normally distributed and required no iterative procedure to be generated. Note that 
by choosing to work with gauge-invariant variables we ensured automatic fulfillment of 
Gauss' law by initial configurations. Having generated an initial configuration, we per- 
formed transformation to a gauge-dependent set of variables as in (§) under an initial 
gauge condition arg0j = 0, Ajj^i = —bjj^i/a and integrated the Hamiltonian equa- 
tions of motion in terms of this latter variable set. We used the fourth-order leapfrog 




algorithm [|12[. In case of AHM and other gauge theories leapfrog algorithms offer, be- 
side their superior energy-conserving properties, the extra benefit of exact local charge 
conservation. For the smallest lattice spacing used, a = 1/16, we varied the time step 
from 0.02 at low temperatures to 0.005 at high temperatures, keeping it small enough to 
suppress any time-step dependence of observables. As a rule, we set the physical size Na 
of the system to 50, known to be large enough to avoid noticeable finite-size effects on 
the rate and other measured observables ||^ . Finally, the value of ^ was 10 throughout 
our simulations. 
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